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DEPTH-RESOLVED SPECTROSCOPIC 
OPTICAL COHERENCE TOMOGRAPHY 

CROSS REFERENCE TO RELATED 
APPLICATIONS 

This application claims priority under 35 U.S.C. 119 from 
Provisional Application Ser. No. 60/048.237. filed Jun. 2. 
1997. the entire disclosure of which is incorporated herein 
by reference. 

BACKGROUND 

Optical imaging of a biological specimen has always been 
a formidable and challenging task because the complex 
microscopic structure of tissues causes strong scattering of 
the incident radiation. Strong scattering in tissue at optical 
wavelengths is due to particulate scattering from cellular 
organelles and other microscopic particles, as well as to 
refractive index variations arising within and between cell 
and tissue layers. 

For over a century, conclusive diagnosis of many diseases 
of cellular origin (such as cancer) has been performed by the 
process of excisional biopsy, comprising the identification, 
removal, histological preparation and optical microscopic 
examination of suspect tissue samples. Many developments 
have taken place to aid the pathologist in interpretation of 
histological microslructure, primarily the development of a 
wide variety of histochemical stains specific to the biochem- 
istry of tissue rnicTostructures. This technique provides 
sufficient resolution to visualize individual cells within the 
framework of the surrounding gross tissue structure. In the 
last several years, a revolution has been stimulated in the 
field of ultra-high resolution microscopy for biomedical 
applications. Ultra-high resolution microscopy allows visu- 
alization of sub-cellular and sub-nuclear structures. This has 
resulted in the invention of tools for high resolution optical 
im a gin g, including near field scanning microscopy, standing 
wave fluorescence microscopy, and digital deconvolution 
microscopy. 

These technologies are primarily designed for imaging 
features at or near the surface of materials. Inhomogeneties 
of the refractive indices inside a biological specimen leading 
to multiple scattering limit the probing depth of these 
techniques. Thus, considerable effort is required to cut and 
preserve the samples in order to prepare the specimen to the 
requirements of the microscope. In medical applications, 
this means that suspect tissue sites identified using mini- 
mally invasive diagnostic technologies such as endoscopy 
must still be acquired and processed via routine histological 
examination. This step introduces significant delay and 
expense. 

The invention of confocal microscopy and its advanced 
development in the past few years have provided the 
researcher the capability to study biological specimens 
including living organisms without the need for tissue 
resection and histological processing. However, the pres- 
ence of multiple scattering in samples limits confocal 
microscopy to specimens which are thin and mostly trans- 
parent There is a need, therefore, for new optical methods 
capable of in vivo imaging deeper inside highly scattering 
tissues and other biological specimens. 

Optical coherence tomography ("OCT") is a technology 
that allows for noninvasive, cross-sectional optical imaging 
in biological media with high spatial resolution and high 
sensitivity. OCT is an extension of low coherence or white- 
light interf erometry. in which a low temporal coherence light 
source is utilized to obtain precise localization of reflections 



2.480 

2 

internal to a probed structure along an optic axis (i.e.. as a 
function of depth into the sample). In OCT. this technique is 
extended to enable scanning of the probe beam in a direction 
perpendicular to the optical axis, building up a two- 
5 dimensional reflectivity data set. used to create a cross- 
sectional, gray-scale or false-color image of internal tissue 
backscatter. 

Many studies have suggested the use of elastic backscatter 
or reflectance as a noninvasive diagnostic tool for early 
io detection of several human diseases, including cancer. The 
use of backscattered light is based on the fact that many 
tissue pathologies are accompanied by architectural changes 
at the cellular and sub-cellular level, for example the 
increase in the nuclear to cytoplasmic volume ratio accom- 
15 panying neoplastic conversion. In the near infra-red (NIR) 
zone, the elastic scattering properties of the tissue are most 
strongly affected by changes in tissue features whose dimen- 
sions are on the same order as the NIR wavelength. Pre- 
liminary success in diagnosing cancer in the bladder, skin, 
20 and gastrointestinal tissues has been reported with tech- 
niques based on elastic backscatter spectroscopy. However, 
currently implemented spectroscopic systems do not incor- 
porate depth resolution and thus cannot provide information 
on the degree of infiltration or cancer staging. Although 
25 elastic backscatter spectra can be collected with confocal 
techniques, the turbidity of biological samples in combina- 
tion with the point spread function of confocal microscopes 
limit the penetration depth for acquiring spatially selective 
spectra to no greater than a few hundred micrometers in 
30 most tissues. Many tissue samples have features of interest 
located at a depth more than mat can be probed by confocal 
techniques, but less than that of other sub-surface imaging 
modalities such as ultrasound. Accordingly, there is a need 
for a spectroscopy system mat is capable of obtaining 
35 depth-resolved clastic backscatter spectra from a sample. 
Inelastic scattering processes including fluorescence and 
Raman spectroscopy have also been exploited for noninva- 
sive disease diagnosis. Unlike elastic scattering events, in 
which the incident and scattered radiation are at the same 
40 frequency, in inelastic scattering events all or part of the 
incident optical energy is temporarily absorbed by the atoms 
and/or molecules of the subject tissue, before being remitted 
at a different (usually lower) optical frequency. Thus, inelas- 
tic scattering processes serve as mtimate probes of tissue 
45 biochemistry. Several studies have reported on laser-induced 
fluorescence spectroscopy as a potential early cancer diag- 
nostic in the skin, breast, respiratory, gastrointestinal, and 
urogential tracts. Additional studies have reported on the 
more biochemically specific Raman spectroscopy for char- 
so acterization of atherosclerotic lesions in the coronary 
arteries, as well as for early cancer detection in the gas- 
trointestinal tract and cervix. In all studies of fluorescence, 
Raman, and other inelastic scattering spectroscopies in 
human tissues to- date, means have not been available to 
55 resolve the depth of the scattered signal with micrometer- 
scale resolution. There is thus a need for a spectroscopy 
system capable of obtaining depth-resolved inelastic back- 
scatter spectra as well as elastic backscatter spectra from a 
sample, such as could be obtained by extending Optical 
60 Coherence Tomography to detect inelasticalfy scattered 
light The depth resolved elastic and inelastic backscattering 
spectroscopic information could aid in the detection of the 
shapes and sizes of lesions in an affected organ and could 
thus assist in accurate staging of diseases such as cancer. 
65 The inelastic scattering spectroscopies based on sponta- 
neous fluorescence and spontaneous Raman scattering 
which have been used in medical diagnostic applications to 
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date are not suitable for combination with Optical Coher- the sample has been altered for each acquisition (eg., by 
ence Tomography because they are incoherent scattering squeezing or stretching the tissue sample) or where the 
processes, and thus the scattered light would not be detected sample arm is repositioned slightly for each acquisition; (c) 
with OCT. However, coherent inelastic scattering processes obtaining an auto-power spectrum by performing a Fourier 
do exist, in particular the process of stimulated emission is 5 transform on auto-correlation data; (d) obtaining a cross- 
the coherent analog of spontaneous emission, and stimulated power spectrum for the windowed portion of each cross- 
Raman scattering is the coherence analog of spontaneous correlation data by performing a Fourier transform on the 
Raman scattering. Other stimulated coherent scattering pro- windowed portion of each cross-correlation data set: (e) 
cesses also occur which may find future application in obtaining a transfer function from the ratio of each cross- 
medical diagnostics, for example coherent anti-Stokes io power spectrum to the auto-power spectrum; (f) squaring 
Raman scattering (CARS) and four-wave ixiixing (FWM). each transfer function obtained in step (e) and (g) averaging 
All of these coherent inelastic scattering processes require the magnitude of the squared transfer functions to reveal 
the presence of pump energy which is converted into signal backscattering characteristics of scatterers resident within 
energy in a coherent gain process. Thus, by virtue of their that window. 

coherence, stimulated coherent gain processes are suitable is Based upon the model described above, any form of 

for combination with Optical Coherence Tomography to coherent spectroscopy can be performed in a depth-resolved 

allow for depth resolution of the location of the inelastic manner. Accordingly, the present invention also provides a 

scattering events. Therefore, there exists a need for a system system and method for j^orming stimulated-emission 

which allows for this combination. spectroscopic optical coherence tomography (SE-SPOCT). 

20 Such a system or method includes the means for, or steps of 
SUMMARY directing an intense pump laser at the appropriate frequency 
The present invention provides a technique for depth- to induce md frequency-dependent gain in the 
resolved coherent backscatter spectroscopy. This technique sample volume interrogated The depth resolved spectrum 
is an extension of OCT technology, ILS. patent application obtained according to the steps discussed above will thus 
Ser. No. 09/040,128, filed Mar. 17, 1998, the disclosure of 25 contain features corresponding to the frequency-dependent 
which is incorporated herein by reference, describes an round-trip gain experienced by the OCT source radiation 
improved OCT system mat utilizes a transfer function (inelastic backscattering characteristics of the scatterers resi- 
modeL where the impulse response is interpreted as a ^nt within the window). In a detailed embodiment of 
description of actual locations of reflecting and scattering SE-SPOCT, the pump laser is cycled on and off and gated 
sites within the tissue. Estimation of the impulse response 30 detection is performed to separate the dastic backscattering 
provides the true axial complex reflectivity profile of the characteristics from the inelastic backscattering characteris- 
sample with the equivalent of femtosecond ranging resolu- tics * Additionally, the pump laser may be modulated at a 
tion. An interferogram obtained having the sample replaced certain frequency and gated detection is performed to sepa- 
with a mirror, is the auto-correlation function of the source rate the elastic backscattering characteristics from the inelas- 
optical waveform The interferogram obtained with the 35 tic backscattering characteristics, 
tissue in the sample is the measured output of the system and Similar to SE-SPOCT, the present invention also provides 
is known as the cross-correlation function. By deconvolving a system or method for performing stimulated Raman scat- 
an impulse response profile from the output mterferometric tering spectroscopic OCT (SRS-SPOCT). In this 
. signal, a more accurate description of the tissue sample is embodiment, the system or method includes the means for, 
obtained 40 or the step of directing a high intensity pump light into the 
This model assists in calculating the spectral characteris- sam P Ic infraction region. The depth resolved spectrum 
tics of optical elements over the bandwidth of the source by obtained according to the steps discussed above will thus 
analysis of the interference resulting from internal tissue contmn locaUzed I*** 5 m * e <kpth-resolved backscatter 
reflections. This model may also be extended to spectrally 45 s P cctlum F™** ^stantial vftrationaiyrotational 
analyze light backscattered from particles in turbid media. s P ectral formation of the scatterers within the sample. 
Because individual scatterers in a turbid specimen may be Accordingly, it is an object of the present invention to 
considered as being essentially randomly distributed in provide a system and method for acquiring depth-resolved 
space, the ensemble average of transfer functions obtained backscatter spectra of a sample, utilizing OCT It is a further 
from cross-correlation data windowed to a specific region ^ object of the present invention to provide a system and 
within the sample reveals the backscattering characteristics method for performing stimulated-emission spectroscopic 
of the scatterers localized to that region. From this modeL it OCT. And it is a forther objective of die present invention to 
is determined that the squared magnitude of the frequency provide a system or method for performing stimulated 
domain transfer function correlates with the backscatter Raman scattering spectroscopic OCT. These and other 
spectrum of scatterers. objects and advantages of the present invention will be 
Accordingly, the present invention provides a system or W""* * om f f J** description, the attached draw- 
method for determining depth resolved backscatter charac- m ^ s and <■» claims, 
(eristics of scatterers within a sample which includes the BRIEF DESCRIPTION OF THE DRAWINGS 
means for, or step of averaging (in the Fourier domain) the 

interferogram data obtained over a region of the sample. 60 1 is a block-diagram representation of a conven- 

In one embodiment of the present invention, the system or tional Micheison interferometer; 

method includes the means for, or steps of: (a) acquiring FIG. 2 is an interferogram of a cover slip placed on a piece 

auto-correlation interferogram data from the Low- of ***** b V Michelson interferometer of FIG. 1; 

Coherence interferometer, (b) acquiring multiple sets of FIG. 3a is a flow-diagram representing a method of the 

cross-correlation interferogram data from the Low- 65 present invention; 

Coherence interferometer having the sample under analysis FIG. 36 is a flow-diagram representing a method of the 

in the sample arm, where the distribution of scatterers within present invention; 
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FIGS. 4a-b is a flow-diagram representing a method of detected electric field 26. which is directed to the optical 

the present invention; detector 28. Because a low coherence light source 10 is used. 

FIG. 5 is a schematic, block-diagram representation of an an interferometric signal is produced at the detector 28 when 

OCT data acquisition and digital signal processing system of the sample probe path distance to a reflecting or scattering 

the present invention; 5 site within the sample 20 matches the reference arm length. 

FIG. 6 illustrates experimental result* of a method of the to within a source coherence length. For every reflecting or 

present invention' scattering site within the sample, a fringe pattern will appear 

- * . ^ in the interferometric signal similar to that as shown in FIG. 

FIG. 7 is a schematic flow-diagram representing an alter- 2 . The axial profiles of backscatter versus depth are mea- 

nate embodiment for obtaining a cross-power spectrum lQ surc d by translating the reference mirror 24 and by synchro- 

accordmg to the present invention; flously rec0 rding the envelope of the mterferometric signal 

FIG. 8 is a schematic flow-diagram representing another a t the detector 28. This profile is known as the OCT A-scan 

alternate embodiment for obtaining a cross-power spectrum c f the sample. Two dimensional cross-sectional imaging of 

according to the present invention; the sample is performed by laterally scanning the sample 

FIG. 9 is a schematic flow-diagram representing yet 15 probe 18 during successive A-scans. The resulting data set is 

another alternate embodiment for obtaining a cross-power processed in a computer 30 and displayed as a gray scale or 

spectrum according to the present invention; false color image. A series of two-dimensional images can 

FIG. 10a is a flow-diagram representing a method for be acquired by scannmg me probe beam peq>endicular to the 

obtaining a spectrum of light in the sample arm according to direction of lateral scanning. The series of two-dimensional 

the present invention- 20 images can then be rendered into a three dimensional display 

FIG. 106 is a flow-diagram representation of an alternate or a V xado teee ^^ional display in gray scale or fake 

method for obtaining a spectrum of light in the sample arm C0 L?" . t .„ . t 

according to the present invention; Those of ordinaryskm in the art wiH recognize that 

. . , , . although it is preferred to scan the sample probe 18 with 

FIG. 11a is a schemata block diagram representation of ^ f to ^ le ^ ^ le & ^ scanoed 

an alternate data acquisition and signal processmg system of ^ rcspcct to a stationary sample probe. It also to be 

the present invention; understood, that for the purposes of this disclosure, the term 

FIG. 116 is a schematic block diagram representation of "optical" is to pertain to all ranges of dectromagnetic 

an alternate data acquisition and signal processing system of radiation, and preferably pertains to the range of 100 nanom- 

the present invention; 30 eters to 30,000 nanometers. 

FIG. 11c illustrates a modification of the systems of FIGS. In developing the present invention, a unique transfer 

11a and lib; function model has been developed for OCT interaction with 

FIGS. 12o-€ is a schematic block diagram representation sample, where the impulse response is interpreted as a 

of an alternate data acquisition and signal processing system description of the actual locations of the reflecting and 

of the present invention, including example diagrams rep- 35 scattering sites within the sample. Based upon this model, 

resenting signals obtained by the system; mc transfer function of the system can be calculated from the 

FIG. 13 is an iUustration of the operation of the alternate ^^^TT, s P C< ^ m J nd me <*<f-P™<* *P«*™r 

embodiments of FIGS. 11a, lib. and 12a; of * c ^^fields m the reference and sample arms The 

„ . . „ % , _ estimation of the impulse response from the transfer func- 

HG. 14a is a schematic : block magram representation of m tion mc ^ ^al complex reflectivity profile of me 

a system for gated <x^^s^^ f <* ™ m sample with the equivalent of femtosecond temporal reso- 

me embodiments of FIGS. 11a, 11*, and 12a; tation vxte^ ^ jnodeL it is determined mat the 

FIG. 14b is timing diagram representing the operation of squared magnitude of the frequency domain transfer tunc- 

the system of FIG. 14a; tion correlates with the backscatter spectrum of scatterers 

FIG. 15a is a schematic block diagram representation of 45 within the sample. In particular, because individual scatter- 

another system for gated or synchronous detection for use ers in a turbid specimen are randomly distributed in space, 

with the enibodiments of FIGS. 11a, 11&, and 12a; and the ensemble average of transfer functions obtained from 

FIG. 15b is a timing diagram representing the operation of cross-correlation data windowed to a specific region within 

the system of FIG. 15a. *h e sample reveals the backscattering characteristics of the 

50 scatterers localized to that region. 

DETAILED DESCRIPTION H Model for Low Coherence Interferometry in Thick Scat- 

L Michelson Interferometer taring Media 

As shown in FIG. 1, a conventional scanning Michelson The present invention is based on a systems theory model 

interferometer can be utilized to obtain the depth resolved which treats the interaction of the low coherence interfer- 

rneasurements of reflectors and scatterers in a sample. A low 55 ometer with the specimen as a linear shift invariant (LSI) 

coherence light source 10 is separated into two beams by a system. For the scan lengths of a few mm typical of OCT 

50/50 beam splitter 16, fifty percent of the light power is imaging, it is assumed that group velocity dispersion is 

transmitted to a sample arm 12 and the rernaining fifty negligible and the group and phase velocities are equal. In 

percent is directed to a reference arm 14. The sample arm 12 the NIR region absorption causes negligible attenuation as 

includes a sample probe 18, which focuses the sample beam 60 compared to the attenuation due to multiple scattering within 

into the sample 20 and collects the retrorcflected light from tissues. This model does not take into account the attenua- 

the sample. The reference arm 14 includes a reference probe tion of light due to multiple scattering as well as absorption. 

22 which transmits the reference beam onto a retroreflecting An optical wave with an electric field with space and time 

rnirror 24, translating towards or away from the reference dependence expressed in scalar form as 2e t (t z) is assumed 

probe, and collects the light retrorcflected back from the 65 to be incident on a Michelson interferometer, as illustrated 

mirror 24.The rerroreflected beams from the sample 20 and in FIG- 1- An arbitrary scaling factor of 2 has been intro- 

mirror 24 are combined again in the beam splitter 16 into a duced for the sake of convenience. The electric field can be 
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expressed as 2e,(t. z^e^t. z)cos(2 Jtf 0 t-koZ+«Kt. z)) where Note that R„<A1). and R,,<A1) are the complex analytic 

to is the central frequency of the source spectrum and representations of the signals measured, i.e.. only the real 

koWloJ ' is the corresponding wave number. ]<, is the ^ of r <A1) and £ (A1 actuaU 

Z Th "! 8 T^f - HCre ' Zl °^ * C We "I™** Held interaction with the LSI 

time and distance traveled by the wave, respectively. 2e Q ,(t. 5 , , . • 

z) and «zj) are the time and space dependenTamplitude and aS 3 ****** Unction H(k) (where k is wavenumber) 

phase of the wave. The analytic signal representation of this J 1 ** inVCfSe FoUriCr transform is bn P ulse res P° nse 

electric field is given by 2e 0 ,(t, z)exp[j(2 K^t-^^t, z))] h(z) ' 

where j=^(-l). In most cases we will use this complex «X-2>=?x-c)®ff*(c) (b> 

analytic representation to denote even those quantities that io 

rZ^ f °L $imPliCit ?: T aXky 2eXt Z)=2CoXt Z)CXP wherc • represents convolution. Note that shift invariance 

D«f.z)]is Wn as the con^lex envelop of the ^ omission of mc tcrms ^ m ^ jon. ^ 

can wnte the term 2pt^ t-k* z as k^vt-z) where v is the convolution theorem leads to 
group velocity at the source center frequency. Here we 

assume that the group and phase velocities are equal. We 15 £,(*)=£,(*)#*(*) (c) 
assume that the complex envelope is varying much slowly 

on the scale of the source wavelength. We further require the where E f (k). E/k), and H(k) are the Fourier transforms of 

T^u^T fJ* f fo ° Cti °? 0f (vt ~ z) ' 111135 11 is e/z),eXz).andh(z),res^^ 

denoted by 2e,<vt~z)=2e 0l <vt-z)expW>(vt~z)]. All the tern- function. The LSI assumption provides: 

poral dependence is expressed in terms of a path delay and 20 97 tt 

thus electric fields can be treated as pure functions of space. e/w-^/^,)^^)) (d) 

Apart from the center frequency, all the information about Inserting Eq. (d) in Eq. (a) leads to 

the wave can be derived from mis complex envelope. 

Therefore we will use the term electric field to denote the ^mHUM^U «^»» (e) 

field e itself as well as its complex envelope e. In the case 25 ^(ky^jktf^k) ^ 

when the group and phase velocities are unequal. ^ 

c(vt-z)=2e 0 Xvt-z)ex^ where c is Note that according to the Wiener-KhiiK&in theorem 

the phase velocity. This analysis will still be applicable. At Fourier transforming the autocorrelation and cross- 

the 50/50 beamsplitter, half the power gets coupled into the correlation functions gives us SJk) and S 6 <k), respectively 

reference arm and the rest into the sample arm. Since the *> s*(k) and S*<k) are the auto-power and cross-power spectral 

power is proportional to the square of the electric field, the deities, relectively. 

electric field at the beam splitter in the reference arm is given This analysis will be useful when raw interferograms are 

by V2e t <vtX setting zM) at the beamspHtter. THe lightwave mcasorcd £ h ^vcrnc* to measure the com^nve- 

readung the reference mirror assumes the form V2e{vt-4A , * - • . . ^ ^ 

where i; and 1, are the optical path lengths in the reference * lopes of the interf erometric data by demodulating R,<Al), 

and sample arms, respectively. If the mirror is ideal it will Sucfa a demodulation assists in noise reduction 

reflect back the same electric field. Then the field returning which °^«wise needs to be achieved by using a bandpass 

to the beam splitter is indicated as V2e{vt-21 ) >/2e<vt-l ) witn sharp cutoffs. Demodulation also allows use of 

is the field reaching the sample. In the sample arm, the light iower sampling frequency while digitizing the data. R^Dl) 

interacts with the specimen and the backscattered field is 40 and RJJ>T) are complex envelopes of R^Al) and R^M), 

given by vTe/vt-l,). The backscattered field reaching the respectively. Thus R^l) and 1^(01) win be measured after 

beam splitter takes the form V2e/yt-21 s ). coherent demodulation. Note that tR*(Dl)| and \RJPI% will 

The fields returning from the reference and sample arms °e nieasured after incoherent demodulation, 

again get separated 50/50 into the arms consisting of the R^(D1) is obtained by coherently demodulating R^(A1) at 

source and detector. Therefore the interference of e, (vt-2^) 45 the center wavenumber ko. Note that when implernented in 

and e, (vt-2i,) is incident on the detector. The reference arm hardware, the demodulation frequency needs to be specified 

length can be varied by various means. One such method is m tcrms of the temporal frequency of the detector current 

mechanically scanning the reference mirror causing sweep- The center frequency of the detector current is 
ing of the reference arm length 1^ Since the detector 

response time (e.g., nanosecond to microsecond for typical 50 *o 

optical receivers) is much longer than die optical wave /r ~ v *2* 
period (-10 15 second) the photocurrent generated by a 
square law detector is proportional to 

which is same as a Doppler shift frequency on the reference 

ts-R^M^MeXYt+Mp <ai) 55 arm light is me scan rate of optical phase deky. V^V r 

The superscript * indicates a complex conjugate. An inter- ror a simple mechanical reference mirror translator. The 

ferogram obtained having the sample replaced with a mirror demodulation frequency can be chosen to be equal to f r for 

in the sample arm, is the "auto-correlation function" R^Al). demodulating the detector current to obtain R^Dl) and 

The tissue is modeled as an LSI system. The "auto- fsX 1 ^)- 

correlation function" R*<A1) of the source optical waveform 60 , Tte f ? e lc f s aQai J^ me s /^ m CX3Sa ^ tx env< ^ 

is treated as an input to the LSI system The interferogram lopcs of * c clect " < l fi f lds md rcs P° nsc ' We 

obtained with the sample in the sample arm is the measured * C T"* 1 ficid cnvclo P e faction with 

output of the LSI system, known as the "cross-correlation ^ san f C f * * msff * wncti . 011 H *> < whcrc k * 

fil „ 5 /An ™ _ A , . . ^ wavenumber) whose inverse Fourier transform is the 

function R^Al). The cross- and auto^oirelation functions 65 impulse response h(z): 

are expressed as: ' 

RAM)^i#>e,*(y*<*rp, k^y^elvf^^t!^ (a2) V2e/ Z )=*2cfc)®h*(-<;) and V2£ # (ik)=V2£ < (*)ff«(Jt) (U) 
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p/JT 5 ® P ? reS T S C COn 7 ,U f? T^T" ^ T* wherc a Fourier fransfo ™ * elated by A . 

J(k) are Founer transforms of e,(z) and e/z). respectively. to practice . me correlatio £ functions ^ in ^ aJ a2 

Note that hCvt-zHtfYt-zjexpUd^vt-z)] and and 3 are hard to measure. The measured (or estimated) 
HCk^Hflc+ko). Also the LSI assumption leads to correlation functions are influenced by the properties of the 

>/2ej(vi-t,y=^2f^i-i,)9h(-(vt-i,)) (ib) 5 optical elements, the measurement electronics, and data 

acquisition systems, and various noise sources. Therefore 

„ , . . „ . £ . » what we measure are "estimates" of R„(Al)(orR„<Al)) and 

The fields returning from the reference and sample arms £ , A i Vor R / Am Hnu ,. v „ fnr ,„ . . .. " '. . 

again get separated 50/50 into the arms consSufg of the ^^SS^ 

source and dcteaor. Therefore the interference of e, (vt-2i r ) 10 the symbols g^Xcr R <A1)) and ^ 

and e (vt-21,) is incident on the detector. The reference arm indicate the "estimates" of auto-correlation and ^ross 

length can be vaned by various means. One such method is correlation functions, respectively. For the purposes of 

mechanically scanning the reference mirror causing sweep- clarity, the terms "autocorrelation and cross-correlation 

rng of the reference arm length \. Since the detector functions" may be used for the terms, "the estimates of 

response time (e.g., nanosecond to microsecond for typical 15 auto-correlation and cross-correlation functions" in describ- 

optical receivers) is much longer than the optical wave ing the present inventions. 

period (~10~ 15 second) the complex envelope of the photo- Similarly the measured (or estimated) power spectra are 

current generated by a square law detector is proportional to influenced by the properties of the optical elements, the 

measurement electronics, and data acquisition systems, and 

'D^{e/i*-yr>^A*-XMe^v*}+< f vt-Zi,)}*>, =*l<£vtY-* f (vt+ 20 various noise sources. Therefore what we measure are 

A)lW*>^vW>, (2) -estimates" of S,,<k) (or S,<k)) and S u (k) (or S*(k)). 

where o denotes integrating over the detector response ^ forthe of ^^scopy algorithms and 

time which is long compared to me electric field period, and < ? aims we wlU sdU usc me ^bols SJti (or S*<k)) and 

2(l,-l r )=D i is the round trip optical path length difference, S «® ( or S*(k)), to indicate the "estimates" of auto-power 

When the reference arm length is scanned at a constant 25 cross-power spectra, respectively. For the purposes of 

velocity, after filtering out the dc components, the time clarity, the terms "auto-power and cross-power spectra" may 

varying components of Eq. 2 reduce to ^ uscd mc toms, "the estimates of auto-power and 

cross-power spectra" in describing the present inventions. 

^~*^ao^"K(+a/> (3) While we describe the specific case of a device which u ses 

, . , , . 30 infra-red light source, the spectroscopy procedure is appli- 

which is just the cross-correlation between the complex cable to any mterferometric device muminated by any 

envelopes of the fields returning from the reference and electromagnetic radiation source. 

sample arms, i n is the complex envelope of the correspond- Eqs. b, c, la, lb, and 5 we describe the light-specimen 

ing current at the photoreceiver output We can also obtain interaction as a linear shift invariant system We describe the 

tirc autocorrelation function of the source field Rjpi) by 35 deconvolution methods based on Eqs. e£ and 6. It should be 

performing the same operation with a imnor in me sample apparent to a person skilled in the art that the interaction 

Jj£ °T h(Z)= ^ Z) t^M^d^C; described by Eqs, b, c ei, la, lb, 5, and 6 can be exploited 

(vt)e,*(vtHAl)>. For an optical source with a well charac- by many other methods in space/time domain as well as 

tenzed spectrum, the form of the autocorrelation function frequency domain including iterative deconvolution 

R^l) is calculated explicitly by computing the inverse 40 methods, etc. This model also forms the basis of "blind" 

Fourier transform of the power spectral density. For a deconvolution methods which do not use a priori informa- 

superlurmnescent diode source ar^oximated by a Gaussian tion about the auto-correlation function but assume mat it 

power spectrum, we obtain convolves with the impulse response. 

The true transfer function H(k) is rarely estimated or 

Real part of J^a/) ~ M^U " t fcJ-co^Dft (4) 45 mcas ?^J* P^^ad cases, what we get is an "esti- 

1 U) J * mate" of the transfer function which is different than the true 

l c = A{]D2)/{Dk) transfer function. One can obtain this estimate in various 

ways. One such method is taking the ratio of cross-power 
spectrum and the auto-power spectrum and taking the com- 
where is the central wavenumber of the source and D k is 50 P lcx conjugate of the ratio. 

the FWHM spectral width. The envelope of the detected Next we interpret the impulse response h(z) and the 
signal (plotted as a function of (lM r )) from a reflection in transfer function H(k). All information regarding the spa- 
the sample arm is a Gaussian function centered at zero tiallv varying and frequency dependent complex reflectivity 
reference arm delay and with a FWHM width equal to the °* mc sample contained in the space domain function h(z). 
coherence length i 0 . Substituting Eq. 1 in Eq. 3 results in 55 Tne tnterpretation for deconvolution is an approximation of 
RU^^R^D^hiDi). (5) mat spectroscopy. In the case of deconvolution, we 

model tissue as a body having several layers of materials 
possessing different refractive indices. The impulse response 
According to the Wiener-Kmnchin theorem, Fourier trans- can be mterpreted as a description of the actual locations and 
forming the autocorrelation and cross-correlation functions 60 reflectivities of reflecting sites within the sample arising 
givesusS^k)andS^(k).S 6 <k)andS^k)aremeaut<H>ower from index of refraction inhornogeneities. The impulse 
and cross-power spectral densities, respectively. We form an response takes a form of a series of spikes (i.e., delta 
estimate of the transfer function H(k) in an arbitrary turbid functions) which are located at the reflection sites while 
sample which may contain many closely spaced reflections dealing with discrete data. These spikes have heights equal 
by Fourier transforrning both sides of Eq. 5: 65 to the electric field reflectivities at these sites within the 

sjifrsSMHik), Hikyx^kys^k). h(z)*^H(k) (6) tissue. This irrformation is available over a length L within 

the tissue where L is the reference mirror scan length. 
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Inordertodescribetheimpulseresponsequantitatively.it response would be complex because of phase delays 

is convenient to deal with discrete representations of the occurred due to different depths of the particles. Therefore 

impulse response and crosscorrelation function. It will also th e effective H(k) would also be complex. A similar phe- 

be helpful to do so since we measure quantized values of the nomenon would occur if the particles were dissimilar from 

discrete representations of crosscorrelation functions using 5 each ^ ^ ^ ^ (fl(k)|2 wouW flt cffectiye 

a computer. We assume that M samples of R„(D1) are 5ackscattering cross-section of the scattered at a wavenum- 

acquired over a length L. Suppose Az is the sampling . , ,„ « , ~ , l2 . f J , . , , 

interval. The numbers ber k We & oyt tot * H ( k H ls reiated t0 the weighted 

average of the reflectivities of heterogeneous scatterers. If a 

lQ low-coherence light source is used the reflections from only 

**<i), R„(2), . . . , R„(n), . . . , r*(M) tnose scatterers located within a few coherence lengths 

would add coherently. In that case |H(k)p relates to the 
denote M measurements of a cross-correlation function weighted average of the backscatter spectra of these non- 
sampled at distances Az, 2Az, 3Az This implies that homogenous scatterers. 

RiAn^R^n^h/in). (j) J5 Let us examine the complex envelope h(n) quantitatively. 

For a homogeneous medium, we can write 

R jV (n) and h(n) are discrete representations of the autocor- 
relation function and the target impulse response, respec- h(ny=b(n)®c(n) (io> 
tively. Now, as discussed earlier, h(n) assumes a form of a 

series of spikes (Le., discrete delta functions) which are 20 where b(n) is a Bernoulli sequence as discussed earlier and 

located at the reflection sites. In a tissue sample, one usually c(n) is the inverse Fourier transform of C(k). Qk) denotes 

does not find a reflection site at every sample. Therefore the spectrally dependent electric field backscattering coefficient 

probabmtyofoccurrenceofanmterfaceatasampleismuch of the scatterers. We call c(n) as the specific impulse 

less than one. A random sequence of zeros and ones is response of a single scatterer. Thus the frequency dependent 

known as a Bernoulli sequence. If the adjacent elements in 25 backscattering cross-section of particles is represented by 

such a sequence are completely unrelated, then such a |Qk)| 2 . Hie Fourier transform of b(n) is given by B(k) and 

sequence is known as a white sequence (since the power E{iB(k)| 2 } represents the expectation value (Le,, the statis- 

spectnim of such a sequence is white). We could possibly tical average) of |B(k)p and is equivalent to the power 

represent the impulse response by a Bernoulli event spectral density of the Bernoulli sequence. For a white 

sequence b(n). The sequence provides a one every time an 30 sequence, this can be assumed to be equal to a constant (in 

interface occurs and a zero in the absence of a reflector at the this case can be assumed to be 1 for simplicity) for all 

sample point However, the amplitude of the reflectivity is wavenumbers. Fourier transforming both sides of Eq. 10 

not constant and fluctuates randomly due to refractive index provides 
inhomogeneities and hence we need to multiply this 

sequence by a Gaussian random number generator g(n). g(n) 35 /?(*)^*)q*> (U) 



is also a white sequence. Therefore 

*(«>=&(»)*(*>). (8) 



Taking magnitude square of both sides gives, 

moHiKQftakf. <i2) 



This essentially means that when at every point in b(n) a one ^ 

occurs (i.e., a reflection site occurs), we turn on our Gaussian Obtaining ensemble averages on both sides yields, 

randomnumber generator and replace the one by the output Etto&UEttokWrfM m , 

of the generator. The value of this random number represents wwwhzmmicwrh (13) 

the reflectivity at that point If the reflectivities are complex. since C(k) is not a randomly varying function, we get 
one can generate a complex Gaussian random number. We 45 

are assuming mat all reflections occur at the interfaces and ^l^^l^lff^lla^i+Ia^F^Jt)! 2 . <i4) 
the effects of point scatterers do not interfere with the 

process of impulse response estimation. Thus averaging several measurements of |H(k)p over a 

In order to perform spatially resolved spectroscopy, a volume cell of interest can provide the backscatter spectrum 

more complicated interpretation of the transfer function is 50 of ^ individual particles localized within the homogeneous 

required. For a monochromatic wave with an amplitude of medium. While computing the spectra of scatterers located 

one and zero initial phase incident on a reflector, the ^cep inside a specimen, one should note that C(k) has the 

reflected field c s is related to the incident field by the information regarding the backscatter spectrum of scatterers 

relation, as well as round-trip spectral filtering due to the intervening 

55 medium. 

i^txp VK*-z)l i^8(k)£r&(k) cxp \jk(vt-z)] (9) Now as indicated in Eq.(c), H(k) is the quantity that 

actually interacts with the electric field. Using 

where H(k) is the reflection coefficient (i.e., backscattering 

coefficient) at a wavenumber k. H(k) is real for many #(*)^(*+*b) and Cxk^k+^x (i4a) 

particles. Then |H(k)| 2 is the intensity reflectivity (Le., #> £{|ff(*)PH<W- (14b) 
backscattering cross-section) at a wavenumber k. If a mono- 
chromatic beam of light is incident on a group of identical |C(k)| 2 is the quantity that actually interacts with the 
particles located at the same depth, the effective H(k) would electric field and is the actual measure of the backscatter 
just get scaled by the number of scatterers due to coherent spectrum. 

addition of the scattered waves. However, if they are situated 65 In an inhomogeneous medium (such as a tissue 

at different depths, then the reflections from all the scatterers specimen), a mixture of various particles within a few source 

would add coherently (i.e., interfere) and the impulse coherence lengths can be described as 
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ft(n)=*, (n)toc t (n H^(n^c 7 (n)+b 3 (ri)^c 3 (n) (15) 



(20) 



where b,<n) {i is a natural number} is a Bernoulli sequence 5 
describing positions of ith type of scatterers having the 

specific impulse response c,{n). Let b,<n) be white processes. Now B,<k) is a Fourier transform of b,<n): 
Suppose these processes are statistically independent of each 

other. Fourier domain representation of Eq. 15 is ~ <2i> 

10 

Hik^B^C^B^C^B^k^k^ (16) 

This implies that 

where B,<k) and C,<k) are Fourier transforms of b f <n) and 

c,<n), respectively. Taking magnitude squares on both sides 15 WH= <22) 
yields, *, 

£ Elb t {n)tczp{-j2xkn\ = Pi £ cxp(-j2xkn\ = p,<>(Jt). 

M (17) *» > 

\Hikf=£mk)\ 2 \c t m 2 + 

20 Therefore the cross terms become 

BtikKii^kyc^i^+BiikyCiik^ik^ikH ... 

u (23) 

where M is the total number of types of scatterers. Taking 
expectation values on bom sides, 25 

considering the fact mat S(k) is ju st a spike of height one at 
2 " , (18) k^m the case of adiscrete Fourier t^ 

EWM j - ^E[mkf\Q(kf\ + terms are not measurable for we can measure the spectra 

*"* only within the bandwidth of the source and hence we can 

EiBtikyCiikWikKZik)} + e [Btikydik^ikyQik)) + ... 30 rewrite Eq. 18 as 



2 4? il <24) 

Since we are dealing with white processes, the terms E *uW)Pl«£*Ka«! 2 J 

{jBiOOp}^ where K, are constants. Now K,*s are propor- 35 
tional to the probabilities p f *s of the occurrence of one in 

each Bernoulli process as long as p/s are much smaller than which is same as the weighted average of the spectra of 

0.5. Bernoulli processes can be considered white only within different scatterers situated within a few coherence lengths, 

the bandwidth of the light source. Since the average of a Sincc me coherence length of the source is short the depth 

Bernoulli sequence is non-zero and equals to p,, there is a 40 ovcr which me spectral information is obtained can be 

spike at zero frequency in the power spectrum and tbcoreti- controlled by limiting the reference arm scan to a region of 

cally this process is not white. However, since our measure- interest within the sample, or alternatively by digitally 

ments are limited to the bandwidth of the source, we model windowing the region of interest from one or more fill length 

these processes as white. When p/s are small, an increase in reference arm scans. Thus the transfer function H(k) 

p i implies an increase in number of scatterers and hence 45 ob^ed over a wide range of path differences contains all 

stronger echoes are obtained Since there are more zeros available information concerning interaction of the tissue 

than ones, all the higher frequencies have equal strength and specimen with the sample arm light including modifications 

the strength increases with an increase in p f . In the power in ***** am P litude P* 185 ©. 

spectrum, the spike at dc does not increase significantly. Again as indicated in Eq.(c), H(k) is the quantity mat 

Therefore K,- increases. When p/s are close to 0.5 or higher, 50 actually interacts with the electric field. Using 

there are many ones occurring and they may occur next to HCk^H^+ko) and C/k)=€,< k+ ko)< 
each other. That would imply the presence of a scatterer at 

every sarrmle, wWchisrareandunlikdy,andinsuchacase, * * ( 24a) 

we should be sampling more often. Also, even if the reflec- E{\R{kf j = £ JQ|Ck*)| 2 1 

tions are stronger, most of the power is in the dc level spike 55 *■> > 
and an increase in p/s simply implies stronger spikes at dc. 

The lower row on the right hand side of Eq. 18 can be 111 a™?** for a monochromatic wave incident on a 
written as reflector, H(k) is the reflection coefficient (i.e., backscatter- 

60 ing coefficient) at a wavenumber k. H(k) is real for many 
" (19) particles. |H(k)| 2 is the intensity reflectivity (i.e., backseat- 

2^ EiBiikWkWitiCtik)] tering cross-section) at a wavenumber k. If a monochromatic 

beam of light is incident on a group of identical particles 
located at the same depth, the effective H(k) would just get 
65 scaled by me number of scatterers due to coherent addition 
Since b,(n) are statistically independent of each other, so do of the scattered waves. However, if they are situated at 
B,<k) and &*{k) and we get different depths, then the reflections from all the scatterers 
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would add coherently and h( z) would be complex because of scattering spectra. Since we are interested in backscattering 
phase delays occurred due to different optical depths of the 

particles. Therefore the effective H(k) would also be corn- spectra, one needs to average §„(k), i.e.. compute E{E,(k) 

plex. Thus |H(k)p would be related to the effective back- E\(k)}. We denote E{E,(k)E* s (k)} by S„(k). i.e.. S,/k)=E{ 

scattering cross-section of scatterers at a wavenumber k. If 5 ^ j 

a low-coherence light source is used, the reflections from S Th us S (k) is given by 

scatterers located within the coherence length would add ** 

coherently. In that case |H(k)p would be related to the ^ (k) = $ k {kjE\\ftikf\ 

effective backscatter spectrum of scatterers. " * II M 

The impulse response h(z) essentially is a convolution of 10 = [\$ u {kf ) fSM 
two functions, viz.. a function b(z) which describes locations 
of scatterers and a function c(z) which is the inverse Fourier 

transform of C(k). C(k) denotes spectrally dependent elec- for a homogeneous medium, 
trie field backscattering coefficient of the scatterers. The 

wavenumber dependent backscattering cross-section of the 15 = 5^*)£||#(jfc)| 2 ) 
particles is represented by |C(k)| 2 . We label c(z) as the 

specific impulse response of scatterers. While dealing with = SiK*)|a*)| 2 
discrete data, h(z) is represented by h(n). For a homoge- 
neous medium, we can write h(n)=b(n)#c(n) where b(n) and 

c(n) are discrete representations of b(z) and c(z). respec- 20 and for a heterogeneous medium, 

tively. We could represent b(n) by a white Bernoulli Sjfy4jfflL m "KAtXktf 

sequence which is a series of randomly occurring zeros and 1 

ones and adjacent elements of the sequence are not related m Depth Resolved Spectroscopy 

to each other. As shown in Eq. 14, above, because averaging Following from the above model, the present invention 

estimates of jB(k)j 2 for a white sequence is assumed to be 25 provides a system and method for determining depth 

equal to a constant (1 for simplicity), averaging estimates of resolved backscatter characteristics of scatterers within a 

|H(k)p (obtained from different locations in a homogeneous sample by averaging (in the Fourier domain) the interfero- 

section of the specimen) provides an estimate of |C(k)| 2 gram data obtained over a region of the sample. In particular, 

which is the elastic backscatter spectrum of scatterers resid- the backscatter characteristics are obtained according to the 
ing within that homogeneous region. The short source 30 steps as illustrated in FIGS. 3a, 36, 4a and 4b. 

coherence length allows to control the depth over which the As shown in FIG. 3a. a first embodiment of a method for 

spectral information needs to be measured. This is achieved determining depth resolved backscatter characteristics of 

by selecting the reference arm scan to an area of interest in scatterers within a sample is illustrated. As indicated in step 

the sample, or alternatively by windowing the area of 31 ^ auto-correlation function R^Al) over a predeter- 

interest from one or more full length reference arm scans. D ^ acquired from an OCT system having an 

While computing the spectra of scatterers located deep optical reflector in the sample arm; and as indicated in step 

medium. Thus the interferometac signal contains depth 40 . . . . . _ . . 

resolved information about both the spatial distributor of As wfll be abated by those of ordwary skfll infte art 

sc^ng centers widunafcsuesampk. as ^eUasabcutthe ff^^Trf^^^T^^i" 
spectralcharacteristics of the individual scatterers. 311 001 Sysfcm - f of Wlthul of * c 

demodulate the mterferometnc data. Ofk) would be esti- * D . ,x~Z _ a - * ^ *T 

t - \-> t~ t-y specimen itself; the auto-correlaUon function can be mod- 
mated usmg such a system. However, |C(k)P or E{|H(k)P } eled ming me Monnation about the radiation source; and 
are preferably calculated as shown in Eq. 14b and Eq. 24a the auto-correlation function can be also calculated using the 
since they represent true interacting spectral characteristics knowledge of the source power spectral density. For 
of the tissue. If we process interferometricdata directly, then ^ me inverse Fourier transform of the measured 
|C(k)| 2 or E{jH(k)] 2 } will be obtained by performing the source power spectrum would provide an estimate of the 
above analysis. auto-correlation function. Additionally, since the auto-power 
The above model also allows one to measure actual spectrum is nothing but the source power spectrum, the 
spectrum of the light in the sample arm. The spectrum of auto-power spectrum can be obtained using the knowledge 
light in the sample arm is defined as 55 of the source. For instance, the source power spectrum 
§ rttF jvb*jv\ measured using any spectrometer or a spectrum analyzer 
SJ.*>£A ]E A ) would provide an estimate of the auto-power spectrum. 

Using E (k>=E{k)H*(k) we cet M m ^ catcA m ste P s ^ 35 and 36. multiple sets of 

* * ' cross-correlation data are obtained, where the distribution of 

_ 6o scatterers within the sample has been altered for each 

$ = 5M*)|i9<*)| acquisition or where the sample arm is repositioned sligjitly 

- 15 kf /$ ik) for aot 3 u * s * tioiL ^ indicated in step 33 the cross- 

U U correlation function ikj^hX) over a full length L is acquired 

from the OCT system having the biological tissue sample in 
Thus the spectrum of light in the sample arm is the source 65 the sample arm; next, as indicated in step 34. the cross- 
power spectrum multiplied by a function which has infor- correlation function is segmented into short segments of 
matron regarding scatterers 1 locations and scatterers* back- length D each to obtain depth resolved spectra with axial 
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resolution D. Length D is typically longer than or equal to slightly for each acquisition. As indicated in step 33a the 

the source coherence length. These segments of cross- cross-correlation function over a full length L is acquired 

correlation data R'„(A1) are called R„(A1). As indicated in from the OCT system having the biological tissue sample in 

step 35. it is deteimined whether a sufficient number of the sample arm; next, as indicated in step 42 the cross- 

cross-correlauon functions have been acquired; and as indi- s correlation function is demodulated at the center wavenum- 

cated in step 36. if more data ,s needed, the distribution of ber of the source; next, as indicated in step 34a the 

The number of sets of «mJLA. .obtain X"3K£3E Ta^T^ 
depends ;upon the desired estimation accuracy of E{|H(k)| 2 }. 10 is determined whether a sufficient number of cross- 
ed typically ranges from 50 to 500 scans. In step 36, the correlation functions have been acquired; and as indicated in 
distribution of scatterers can be altered, for example, by step 36, if more data is needed, the distribution of scatterers 
lightly squeezing or stretching the tissue sample, by propa- within the sample is altered or the sample arm is sliehtlv 
gating sound or ultrasonic waves through the tissue while repositioned prior to returning to step 33 
acquiring the data, or by repositioning the sample arm with 15 As indicated in step 37a, the cross-power spectrum is 
respect to the sample For the purposes of the present obtained from each of the demodulated cross-correlation 
mvenuon such steps will be referred to as altering the functions by performing a Fourier transform on the cross- 
distribution of scatterers with respect to the sample arm. correlation function. As indicated in step 39a, a transfer 

Jf , in step 36, the sample arm is repositioned, the extent function is determined for each cross-power spectrum by 

of such repositioning depends upon the beam diameter at the 20 taking the ratio of each cross-power spectrum versus the 

depth of interest, and is typically 20-50 mm. The separation auto-power spectrum. As indicated in step 43, each transfer 

between two adjacent scans should be at least the source function calculated is adjusted to shift the axis from the 

beam diameter in the tissue. Since the beam waist is different center wavenumber of the source, thereby removing the 

at different depths within the tissue, the mmimum required effect of the demodulation steps (see Equation 14(a)) And 

separation wfll be different for each depth. 25 as indicated in step 40a, the magnitude of the backscatler 

As indicated in step 37, the cross-power spectrum SJk) spectrum C(k) is calculated by averaging the magnitudes of 

is obtained from each cross-correlation function by perform- the squared transfer functions and then taking the square 

ing a Fourier transform on the cross-correlation function root of the average, 

Ri/Al). As indicated in step 39, a transfer function H(k) is illustrates a slightly modified version of the 

determined for each cross-power spectrum S (k) by taking 30 mcthod illustrated in FIG. 4a. In particular, steps 33a and 

the ratio of each cross-power spectrum § jk) versus the ?£ ™ ****** f? m mG ' <*• As 

a nt^n,.r 01 jtla * j " . ~! mdicated in step 33' the cross-correlation function over a full 

™^r^T^ if ^ m *^ dmStep4(K * e length Lis acquired; however, the full length L is equal to 

magnitude of the backscatter spectrum C(k) is calculated by the desired segment length D and merefore a tferTof 

ZTT*^^*** <> f ^^transfer functions 35 segmenting JSLSS^ Zi t^Z^ 

mG^l^rjZf ^ 2** - , „ As shown in FIG, 5, an OCT data acquisiS system 44 

t^Z^^C^ ^ for performing the above method incluSs a low^herence 

cSfi£^ 33 r d3 1 ^terferometer 46 and, preferably, a cafibration interferom- 

of FIG. 3a are replaced with step 33 in FIG. 3b. As mdicated eter 48. The low-coherence interferometer includes a radia- 

in step 33 1 the cross-correlation function R^(A1) over a full 40 tion source 50, such as a super-lurninescent diode CSLD") 

length L is acquired; however, the full length L is equal to source, a fiber-optic source line 52 coupled between the SLD 

the desired segment length D, and therefore, a step of 50 and a fiber-optic beam splitter (such as a 50/50 fiber 

segmenting the cross-correlation data is not required. coupler) 54. The beam splitter separates the fight received 

FIG. 43 illustrates another embodiment of a method for from the source line 52 into two beams; one transmitted to 

determimng depth resolved backscatter characteristics of 45 a sample arm 56 via an optical fiber 58, and the other to a 

scatterers within a sample. The primary difference between reference arm 60 via an optical fiber 62. The fiber 58 is 

the method illustrated in FIG. 3a and the method illustrated coupled to a sample probe 64 adapted to focus light to a 

in FIG. 4a that the method illustrated in FIG. 4a includes the sample 66 and to receive the light reflected back from the 

steps of coherently demodulating the auto- and cross- sample 66. The reflected light received backfrom the sample 

correlation functions at the center wavenumber of the source so is transmitted back to the beam splitter 54 via the fiber 58. 

prior to obtaining the auto- and cross-power spectrums. Preferably, the sample probe 64 has an adjustable focal 

respectively. Additionally, each transfer function calculated length, thus allowing the adjustment of the focal spot size, 

is adjusted to shift the axis from the center wavenumber of working distance, and depth of focus, 

the source, thereby ridding the effect of the demodulation The fiber 62 is coupled to a reference probe 68 adapted to 

stc P s * . 55 focus the light received from the fiber 62 to a translating 

As indicated m step 31a, the auto-correlation function reference mirror 70 (usually mounted on a galvanometer), 

over a predetermined depth D is acquired from an OCT and to receive the light reflected back from the reference 

system having an optical reflector in the sample arm; as rnirror 70. The reflected light received back from the refer- 

lndicated in step 41 the auto-correlation function is demodu- ence mirror is transmitted back to the beam splitter 54 via 

lated at the center wavenumber of the source; and as 60 the fiber 62. The reflected light received by the beam splitter 

mdicated in step 32a. the auto-power spectrum is obtained 54. back from both the fiber 58 and fiber 62, is combined and 

from the demodulated auto-correlation data by performing a transmitted on the fiber-optic line 72 to the photodetector 74. 

Fourier transform on the demodulated auto-correlation data. The photodetector 74 produces an analog signal 75 rcspon- 

As mdicated in steps 33a, 42, 34a, 35a and 36a, multiple sive to the intensity of the incident electric field An example 

sets of cross-correlation data are obtained, where the distri- 65 of a photodetector for use with the present invention is a 

bution of scatterers within the sample has been altered for Model 1811, commercially available from New Focus, Inc., 

each acquisition or where the sample arm is repositioned Mountain View, Calif. 
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It will be apparent to one of ordinary skill in the art that 
there are many known methods and/or mechanisms for 
injecting the above reference arm delay, other than a trans- 
lating reference mirror. All of these methods, of course, are 
within the scope of the present invention. 

Alternative reference arm optical delay strategies include 
those which modulate the length of the reference arm optical 
fiber by using a piezo-electric fiber stretcher, methods based 
on varying the path length of the reference arm light by 
passing the light through rapidly rotating cubes or other 
rotating optical elements, and methods based on Fourier- 
domain pulse-shaping technology which modulate the group 
delay of the reference arm light by using an angularly 
scanning rnirror to impose a frequency-dependent phase on 
the reference arm light after having been spectral dispersed 
This latter technique, which is the first to have been shown 
capable of modulating the reference arm light fast enough to 
acquire OCT images at video rate, depends upon the fact mat 
the inverse Fourier transform of a phase ramp in the fre- 
quency domain is equal to a group delay in the time domain. 
This latter delay line is also highly dispersive, in that it can 
impose different phase and group delays upon the reference 
arm light For such a dispersive delay line, the OCT inter- 
ferogram fringe spacing depends upon the reference arm 
phase delay, while the position of the interferogram enve- 
lope at any time depends upon the reference arm group 
delay. All types of delay lines can be described as imposing 
a Doppler shift frequency 




on the reference arm light where Doppler shift frequency in 
mis context is defined as the time derivative of the phase of 
the central frequency component present in the interfero- 
metric signal R^Al). is the scan rate of optical phase 
delay in the reference arm. V^is also considered as scan rate 
of the optical phase delay difference between the reference 
and sample arms. This definition of Doppler shift frequency 
encompasses all possible reference arm delay technologies. 
The optical path length 76 of the sample arm 56 is different 
for reflecting and scattering sites at different depths (altering 
the scatterer distribution will change the reflector locations), 
while the optical path length 78 of the reference arm 60 
changes with the translation of the reference mirror 70. 
Recording the detector current while translating the refer- 
ence mirror 70 provides interferogram data, which is the 
optical path length dependent cross-correlation function 
R*(A1) of the light retroreflected from the reference mirror 
70 and the sample 66. Collecting interferogram data for a 
point on the surface of the sample for one reference rnirror 
cycle is referred to as collecting an "A-scan." The A-scan 
data provides a one-dimension profile of scattering infor- 
mation of the sample 66 verses depth. 

The analog interferogram data signal 75 produced by the 
photodetector 74, for each A-scan, is sent through data 
processing scheme 80, designed to perform the steps as 
described above in FIGS, 3a-46. The data processing 
scheme 80 includes an analog-to-digital converter 82 for 
converting the analog interferogram data 75 produced by the 
photodetector 74 into a digital interferogram signal 84. The 
digital interferogram signal 84 is sent to a windowed Fourier 
transform algorithm 86 for obtaining the cross-power spec- 
trum S w (k) data 88. Windowed Fourier transform algorithm 
86 uses a Fourier transform algorithm which is available in 
software libraries in commercially available software pack- 
ages such as Lab VIEW supplied by National Instruments, 
Austin, Tex. 
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The cross-power spectrum data is then sent to a process- 
ing algorithm 90 for calculating the transfer function esti- 
mate H'(k) data 92. The processing algorithm 90 is coupled 
to a memory 94 for storing the auto-power spectrum data 

5 and for storing the multiple transfer functions H (k) used for 
calculating the backscatter spectrum C(k). To obtain the 
auto-power spectrum data, the sample 66 is replaced by a 
mirror 66* and the data received by the photodetector 74 is 
the optical path length dependent auto-correlation function 

io R„(Ai) of the source light generated from the light retrore- 
flected from the reference mirror 70 and the sample mirror 
66\ The analog-to-digital converter 82 converts the analog 
autocorrelation function 75 into a digital signal and the 
Fourier transform algorithm 86 then obtains the auto-power 

15 spectrum S f ,< k) data 88. When the processing algorithm 90 
receives the auto-power spectrum S„<k) data, it stores the 
data in the memory 94. Accordingly, the processing algo- 
rithm 90 will have access to the auto-power spectrum S u {k) 
for calculating each estimate of the transfer function H(k) as 

20 described above. 

Each estimate of the transfer function H(k) 92 is prefer- 
ably obtained by the processing algorithm 90 according to 
the ratio of each cross-power spectrum SJk) versus the 
stored autopower spectrum S 6 {k). 

25 Once all of the transfer function estimates fr(k) 92 have 
been calculated, the transfer function estimates fr(k) are 
transmitted to an algorithm 96 for obtaining the average of 
the squared magnitudes of the transfer function estimates, 

3Q and for calculating the backscatter spectrum estimate |C(k)l 2 
98 by taking the square root of the average as discussed 
above with respect to step 40 of FIGS. 3a and 3b. and step 
40a of FKjS. 4a and 4& An algorithm 96 for use with the 
present invention may be implemented in commercially 
available software packages such as LabVESW supplied by 

35 National Instruments, Austin, Tex. 

Note mat the operations described herein have been 
performed and tested in software using packages such as 
LabVTEW and MATLAB. It is also within the scope of the 
invention mat these operations be performed by using hard- 

40 ware DSP devices and circuitry. For example, the Fourier 
transform algorithm 86, the processing algorithm 90 and the 
algorithm 96 may be performed by hardware devices or 
circuits specially designed to perform the steps as described 
above. Such hardware devices or circuits are conventional 

45 and thus will be apparent to those of ordinary skill in the art 
The backscatter spectrum C(k) 98 is transmitted to a 
computer 100 for comparison against backscatter data from 
a 'normal' tissue that is stored in the database 101. Addi- 
tional applications for the backscatter spectrum C(k) 98 are 

50 given below. The computer 100 also preferably generates the 
control signals 104 for controlling the above process. For 
example, the computer may simultaneously control the 
lateral translation of the sample probe 64 and the translation 
of the reference mirror 70; and the computer 100 may also 

55 provide controls for coordinating the deconvolution scheme 
80. Furthermore, it should be apparent to one of ordinary 
skill in the art, that the computer 100 could contain all or 
portions of the data processing scheme 80, or that the data 
processing scheme could be part of a separate analog or 

60 digital circuit, etc. 

It will be apparent to those of ordinary skill in the art that 
a demodulation step/device may be incorporated to the 
system of FIG, 5 to demodulate the interferogram data 75 or 
84 prior to the processing steps/components discussed 

65 above. For example, see FIG. 4b. 

A preliminary demonstration of depth resolved spectros- 
copy using OCT is provided in FIG. 6. In the figure, the 
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spectrum of the light incident on an interference filter is This spectral information could also be used to comple- 

plotted along with the transfer function (i.e.. spectral pass- m ent the diagnostic capability of OCT. Three dimensionaliy 

band characteristic) corresponding to a double pass through resoive4 spectroscopy could assist in studying the extent of 

the filter. The transfer function is obtained using Eq. 6 of the : „ fi - *V f . . - . . ; - e , 

model developed above from an OCT A-scan ofthe filter. 5 "f" ° D °! a SUCh a f Can ^ * f* 

windowed to the vicinity of the glass-air interface from the S1 8 natures of vanous la y ers of ™nnal and ^normal tissue 

rear side of the filter. Referring to FIG. 6. the spectrum of sam P Ies - *™ pathological states at different layers could be 

light (1) incident on a commercial interference filter (IF) and determined by comparing the spectra acquired at these 

the transfer function characteristic (passband) resulting from Ia y ers witn mose m me database. Two dimensional lateral 

a double pass through an interference filter (2) are plotted. scanning of the tissue samples may provide the degree of 

Both spectra were obtained by separately gating and pro- 10 invasion of diseases accurately. 

Z^^^f^ CT T CtliC *T rCSU i^ g IT* FrCSnel V * UM *> Obtain Spectra at a Given Depth 
reflections at the front and rear surfaces of the filter, respec- 

lively. The measured spectral widths correspond well with _ As Ascribed in section HL the cross-correlation function 

the manufacturer's specifications (SLD: spectral width 47 R'«(Ai) corresponding to a lateral (or angular in endoscopy/ 

nm FWHM; filter: bandwidth 10 nmFWHM single-pass). 15 catheterization applications) position of the sample probe 

IV. Depth Resolved Backscatter Fourier Transform Spec- over a full length L is acquired from the OCT system having 

troscopy Application mc biological dssue & me k and * 

such changes include enlargement and changes in^pes of 20 spectrawim axial resolution D.D is typically 

nuclei in colonic adenoma (which is a precursor to colonic chosen t0 lon S er man <* t0 me source coherence 

cancer) and an increased population of inflammatory cells in length. These segments of cross-correlation data R^(A1) are 

ulcerative colitis. These scatterers in the tissue have dimen- called R^(A1). 

sions comparable to the near infra-red wavelengths ofthe Here we elaborate on a method to obtain segments of 

light source (e.g., 1250 nm to 1350 nm for the SLD in our 25 5. , An - , . A - A . A r Z~T 

laboratory). Therefore the scattering process is best 5^*2 *?T g * ^^fje^ range 

described by the phenomenon of "Mie scattering- for which D "? . sdc< ? C(L M &c l dected startm S de P th ^ the values of 

the backscattering cross-sections of these scattering sites are digitized R'^Al) (ie., R'^(n)) are selected corresponding to 

highly frequency dependent In a typical histopathologic^ depth range D. The selected array corresponding to the depth 

assessment changes in cellular structure are examined 30 range D is called digitized R^Al) (La, R*(n)) and is passed 

These include changes in cellular as well as nuclear sizes into a Fourier transform circuit or algorithm to obtain a 

and shapes and (^storing patterns of the cells or nuclei. power spectrum SJk) for that particular depth range. 

Variations in morphology affect the elastic scattering prop- XT _ , J* . . * ^ F . * 

erties of cells. ThTvarktion in backscattering propel of NcXt * a prcdetermme ^ P° mt starting depth, 

the tissue microstructure could be exploited to diagnose 35 <k 3* h window of RJn) is extracted from R*(n). 

various diseases. Depth resolved elastic backscattering spec- From this next depth window, S u (k) is calculated as 

troscopic information could aid in detecting the shapes and described above. New windows will thereafter be repeatedly 

sizes of the lesions in an affected organ and thus assist in extracted and processed to generate a complete spectral 

accurate staging of diseases such as the cancer, profile for the particular R^(A1). 

Extracellular architecture as well as the sizes and shapes 40 This operation can be sumniarized by the following 
of cells and cellular components exhibit variations in dif- equation (known as windowed Fourier transform (WFT) 
ferent tissue types as well as in different layers within the equation): 
same tissue. The frequency dependence of the elastic back- 
scattering properties of biological materials is closely related f m _ x . 
to these morphological changes. Therefore this spectral 45 S^At^J f fain + mVx}a&~M'mto]w(m*z) 
information can be used in detecting shapes, sizes or refrac- [m^in 
five indices of various particles at different depths in the 
tissue specimen. This can be achieved since the backseat- 

tered spectrum is a function of shape, size and the refractive where w(mDz) is the analysis window through which the 

index of the particles and the refractive index of the sur- 50 sampled space-domain cross-correlation function R'*(n) is 

rounding medium. sriifted, Dz is the sampling interval, and N is the number of 

This spectral information may also provide contrast samples contributing to the local spectral estimate centered 

mechanisms based on differential backscattering spectro- at depth nDz. 

scopic properties of the sites localized within different depth Thus N samples of S^nAz, k) are acquired. S u (nAz, qkf) 

regions of the specimen. Range gated spectra thus obtained 55 , 5 _ A * a \u \i TT tJ > J 

can be used to develop various contrast enhancement 1S measurcd at depth nDz ' Here ^"(NAz) * the 

mechanisms. For instance, an OCT image obtained with s ^^g ^rval in wavenumber domain, and q is an 

complex envelope information can be displayed as a gray mteger numbers 
scale image using amplitude information. It is possible to 

compute spectra resolved at different depths in regular 60 sJnte, -[ N 1*4 sJ«Az, -f * ~ 2 W, 
intervals using complex envelope date. This spectral infor- * \2) T \ \ 2 ) r 
mation can be "color" coded (e.g., one can encode the peaks / (N-i\ ) 
or widths of the spectra in various colors). Thus the gray * ^<*^ $4*^ [ 2 J*' J 
scale image can be supplemented with mis "color" infor- 
mation resulting in a colored OCT image. Better differen- 65 

tiation between the layers in OCT images could be achieved denote N measurements of a the cross-power spectrum 

by using such a color information. sampled at spatial frequencies 
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In our experiments we used a simple rectangular window. 
{ N - 2 \ An alternative implementation of the window is to pad the 

Aft R -| 2 , ... N-point analysis window with zeros on either side, increas- 

ing its length in order to enhance the frequency precision. 

„ , - • 5 Tte size (length) of the window is indicated as NDz which 

Thus the spatial resolution of the spectral estimate is given _ . . u . „ A - ¥ . 

/ , . ,*r^ v . . . ... must be shorter than the entire A-scan length L. While a user 

by the window size (NDz), the larger the window — the . . . , . . _ . . . 

. J it _ . ZL . - may choose the window length as short as he wishes, it may 

lower the spatial resolution. But spectral estimation preci- . J ^ . c A . ^ A - A , 4 ^ 

. j> . - . t i • . . . . r . be apparent to those of ordinary skill in the art that due to the 

sion kr is inversely related to the window size and is given «HF<"=m iu wuac vi uiumoij »m «c «i ui« uu C w wc 

^ interference from the reflectors located withm a coherence 

10 length, the axial resolution is limited by the coherence 

ik r =i/(NAz) length. Therefore it is desirable to choose the window length 

longer than or equal to the coherence length. 

This equation is derived using discrete Fourier transform The above steps for converting the interf erogram data into 

properties. Thus the larger the window-the better the ^ ^ spectnjm ^ ^ ^ perfonned by a bank of 

spectral resolution. narrow-band band-pass filters (NBPF). where each NBPF 

If R fa (Al) is coherently demodulated to obtain the com- passes a par ticular wavenumber along the power spectrum 

plex envelope R ^(Al). then a similar procedure provides wavenumber scale. The outputs of each NBPF may be input 

depth resolved spectra. As described in section KL the directly into the transfer function calculation. This method 

cross-correlation function R'^Al) corresponding to a lateral 20 eliminates the need for the windowed Fourier transform 

(or angular in endoscopy/catheterization applications) posi- ckcuit/algorithm and. may also eliminate the need for the 

tion of me sample probe over a full length L is acquired from coherent demodulation circuit/algorithm. Thus, this method 

the OCT system having the biological tissue sample in the provides faster and cheaper signal processing, 

sample arm and is segmented into short segments of length These NBPF may also be implemented by a bank of 

D each to obtain depth resolved spectra with axial resolution 25 demodulators and a corresponding bank of low-pass filters. 

D, D should be longer man or equal to the source coherence where each demodulator demodulates the data at a particular 

length. These segments of cross-correlation data R^Al) are wavenumber along the power spectrum frequency scale and 

called R'^Al). each corresponding low-pass filter filters the output of the 

Here we elaborate on a method to obtain segments of demodulator to only pass a narrow band of wavenumbers as 
R^(A1). First a starting depth and a window of depth range 30 desircd bv uscrs - 

D is selected. At the selected starting depth, the values of These NBPF can be applied directly to R^Al) or R'«(A1) 

digitized R'^Al) (or R^n)) are selected corresponding to to obtain S'^k) or S^k), respectively. Also, one may 

depth range D. The selected array corresponding to the depth demodulate KJM) at a wavenumber other than ko to 
range D is called digitized R^Al) (or R^(n)) and is passed 35 produce a signal which can be fed to NBPF to eventually 

into a Fourier transform circuit or algorithm to obtain a generate Sf/k) or S^(k). 

Thus in this alternate embodiment of the present 
Next, at a predeternuned point past the starting depth, ^ ncccssity of a windowed Foaricr tn£sform 

another depth window of R^n) is extracted from R^n). ^. to mo ^ acc mc powcr spectra may S (k) or S (k) may 

Rom dm next depth window . S*(k) h calcolated as 40 by J*^, Xk of nLow-band band-pass 

described above. New windows will thereafter be repeatedly fflteIS ("NPBFs^.^tee each NBPF passes a particular 

extracted and processed to generate a complete spectral ^ ^ ^ spectrum fr ^ uency ^ 

profile for me particular RV/A1). As shown in FIG. 7. a bank of NBPFs 130 may be 

^operation can be summarized by the following ^ to mter ^ ra w. (non^lemodulated) photodetec- 

^ * tor signal 75 and produce the cross-power spectrum 88\ 

From the cross-power spectrum 88' the processing algorithm 
c , a v ^ ("v 1 *y fl ^A , , • f . i A J 90 is used to generate a transfer function and the algorithm/ 

x estimate. The bank 130 comprises N NBPFs, each having a 
bandwidth of approximately (V^/NDz). Here is the scan 
S^nAz, qkT) is SJk) measured at depth nDz. rate of the group delay of the reference arm as compared to 
Thus the spatial resolution of the spectral estimate is the sample arm. V^=2V r for a simple mechanically trans- 
given by the window size (NDz), the larger the window — lated reference mirror Such a filter provides spatial fre- 
the lower the spatial resolution. But spectral estimation 55 quency precision given by kf=l/(NAz). The center fire- 
precision kf is inversely related to the window size and is quency of each NBPF corresponds to a particular 
given by wavenumber along the wavenumber scale of a power spec- 

_ trum S^(k), centered at the source center wavenumber ko. If 

*--i/(NAz) is ue cdter frequency of a band-pass filter (BPF), then. 

This equation is derived using discrete Fourier transform 60 me corresponding wavenumber is given by: 
properties. Thus the larger the window — the better the 

spectral resolution. k A x _ v * j + f* ^ 

Several types of analysis windows may be used in the \ v,J v t 

circuit/algorithm, including rectangular. Bartlett 

(triangular). Hamming, Hanning, Blackman windows. It is 65 

well known to those skilled in the art that the choice of For a simple mechanically translated reference mirror, 
window may affect the power spectrum estimation accuracy. V g =Vj=2V r 
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NDz; and the next NBPF below the center NBPF wiJJ be 
f c centered at the frequency f r -f 7 -\yNDz. The wavenumber k 

k = v r n is given by: 

5 / V A f< + ft 

Thus, some of the NBPFs will pass frequencies below the k = U i - + y 2* 

Doppler shift frequency of the reference arm and the rest ^ *' 1 

will pass frequencies above the reference arm Doppler shift 

frequency i r For example, the middle NBPF will be cen- ^ ^ embodiment ±c wifld owed Fourier trans- 

tere K ^ e 'il i CC D X P io ii rCqU f DCy /n t ' «> form circuit/step is not needed. From the cross-power spec- 
next NBPF above the center NBPF will be centered at the _ m A " - niY . - f . 0< > . „ « /[ ftani J f 
frequency f r +V JNDZ; and the next NBPF below the center tn3m f 88 * c ^ocessmg algonthm 90 is used to generate a 
NBPF will be centered at the frequency f r -V JNDZ. In this ^ction and the algonthm/arcuit 96 is used to 

alternate embodiment, the coherent demodulation circuit/ calculate the bacl^catter spectrum estimate As discussed 
device and the windowed Fourier transform circuit/step are above in nGS - *° and »«»f a involution step was 
not needed. 15 performed on the interferogram data, the axis of the transfer 

Alternatively, as shown in FIG. 8, the photodetector signal function must be shifted to remove the effect of demodula- 
75 is demodulated at the center wave number of the source, tion. 

by demodulation step/device 132, to produce analog It will be apparent to one of ordinary skill in the art that 
in-phase *T* 134 and analog quadrature "Q" data 136. A the array of complex NBPFs may be replaced by two arrays 
bank of complex NBPFs 138 are positioned after the coher- 20 of NBPFs, one array for the in-phase data, and one array for 
cnt demodulation step/device 132, where each receives bom the quadrature data; The frequency f,is selected such that the 
the analog in-phase "F 134 and the analog quadrature *Q" NBPFs are centered in the few kHz frequencies. Such 
136 components of R^(Al).The bank 138 comprises NBPFs, NBPFs are cheaper, and more readily available man NBPFs 
each having a bandwidth of approximately V^NDz, where centered in the frequencies near the reference arm Doppler 
the center frequency of each NBPF corresponds to a par- 25 shift frequency i r 

ticular wavenumber along the wavenumber scale of the The banks of NBPF, discussed above, may also be 
power spectrum S^(k) 88" centered at zero wavenumber. replaced by a bank of demodulators and low-pass filters. 
Thus, some of the NBPFs will pass frequencies below zero, where the demodulation frequency is same as the center 
and the rest will pass frequencies above zero. The middle frequency of the corresponding BPF and each corresponding 
NBPF will be centered at zero frequency; the next NBPF 30 low-pass filter has a bandwidth as desired by users, 
above the center NBPF will be centered at the frequency ft should be obvious to one skilled in the art mat the 
V g/NDz; and the next NBPF below the center NBPF will be approach of using a bank of NBPFs or demodulators to 
centered at the frequency— V ^/NDz. In mis alternate obtain spatially-localized frequency information is quite 
embodiment, the windowed Fourier transform circuit/step is general , and is not necessarily limited to the case of using all 
not needed. Rom the cross-power spectrum 88" me pro- 35 NBPFs symmetrically distributed around the reference arm 
cessing algorithm 90 is used to generate a transfer function Doppler shift frequency, or around baseband or at evenly 
and the algorithm/circuit 96 is used to calculate the back- spaced frequencies, or even all with the same pass band- 
scatter spectrum estimate. As discussed above in FIGS. Aa width. The user may have the option to select any frequency 
and 4b, because a deconvolution step was performed on the to monitor, with any bandwidth desired, 
interferogram data, the axis of the transfer function must be 40 VL Method to Measure Actual Spectrum of the Light in the 
shifted to remove the effect of demodulation. Sample Arm 

If i c is the center frequency of a band-pass filter (BPF), Following from the model described above in Section H. 
then, the corresponding wavenumber is given by: the present invention also provides a system/method for 

detennining or measuring the actual spectrum S M (k) of the 
£ 45 light in tiie sample arm 56. As shown in FIG. 10fl, a first 

v ^2* embodiment of this method includes a step 148 of obtaining 

an auto-power spectrum according to any of the procedures 
described herein; a step 150 of obtaining multiple cross- 
It will be apparent to one of ordinary skill in the art that the power spectrum readings from the interferometer signal 
array of complex NBPFs may be replaced by two arrays of 50 according to any of the procedures described herein, where 
NBPFs, one array for the in-phase data, and one array for the the distribution of the scatterers within the sample is altered 
quadrature data; and the square-root of the sum of the w ith respect to the sample arm for each reading; and a step 
squares of the output of each array will yield the resultant | 52 of calculating the spectrum S„(k) of the light in the 
power spectrum S^k). sample arm according to the following equation: 

Finally, as shown in FIG. 9, the interferometer signal 75 55 
can be coherently demodulated at frequency f 7 , as shown in 
block 140, to produce in-phase data 142 and quadrature data 

144. The in-phase and quadrature data is fed into a bank of As shown in FIG. 106. an alternate embodiment of this 
complex NBPFs 146 to generate a cross-power spectrum method includes a step 154 of obtaining a demodulated 
88' The bank comprises N NBPFs, each having a band- 60 auto-power spectrum according to any of the procedures 
width of approximately V/NDz, where the center frequency described herein; a step 156 of obtaining multiple demodu- 
of each NBPF corresponds to a particular wavenumber lated cross-power spectrum readings from the interferometer 
along the wavenumber scale of a power spectrum centered signal according to any of the procedures described herein, 
at frequency f 7 . Thus, some of the NBPFs will pass frequen- where the distribution of the scatterers within the sample is 
cies below f 7 . and the rest will pass frequencies above f 7 . The 65 altered with respect to the sample arm for each reading; a 
middle NBPF will be centered at i,. the next NBPF above the step 158 of adjusting the axis of each of the demodulated 
center NBPF will be centered at the frequency f r -f/fV^/ auto- and cross-power spectra to remove the effect of the 
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demodulation step: and a step 160 of calculating the spec- of the calibration interferogram (eg., zero crossings) using 

trum S„(k) of the light in the sample arm according to the a thresholding or pattern recognition algorithm; and 

following equation: re-sample the SLD interferogram data at the regular inter- 

. vals using interpolation routines. 

SJJl\-E{$M }/S„ (k) 5 vjjj stimulated Coherent Spectroscopic Optical Coherence 

VH. Reference Arm Optical Path-length Calibration Tomogaphy (SC-SPOCT) 

It is advantageous, in the above data processing scheme The present invention also provides a system and method 
that the auto-correiation and cross-correlation functions be for performing three-dimensionally resolved coherent spec- 
measured with the sub-micron accuracy. Therefore to troscopy by taking advantage of the depth-resolving capa- 
enhance the accuracy of the low coherence interferogram 10 OI * OCT, In principle any coherent scattering process 
acquisition, a long coherence length calibration interferom- detected by OCT, as long as the scattered light 
eter 48 is incorporated into the system to accurately monitor remains within the bandwidth of the OCT source light 
and compensate for the inevitable velocity fluctuations of spectrum. In particular, we disclose a method for depth- 
the reference mirror 70. resolving stimulated coherent scattering processes using a 

As shown in FIG. 5, the calibration interferometer 48 15 system based upon the SPOCT concept Stimulated coherent 

includes a long-coherence length, narrow-band laser illumi- scattering processes involve a transfer of energy from pho- 

nation source 106, such as a helium neon (He— Ne) laser or tons * n m intense pump beam and from excited states of an 

a distributed feed-back diode laser (DFB diode laser), a atom °r molecule to a typically weaker probe beam. The 

reference probe 108, and a sample probe 110. The narrow- I* 0 ** photon frequency is typically Stokes shifted (i.e, is at 

band iUumination source must have a coherence length that 20 a iower frequency) from the pump photon frequency, 

is longer than the region (depth) in the sample 66 that is Examples of stimulated coherent scattering processes 

being scanned (for example the He— Ne laser has a coher- include stimulated emission, stimulated Raman scattering, 

ence length of several meters). coherent anti-Stokes Raman scattering (in which case the 

Hie mumination source 106 transmits to a beam splitter P** 0 ** photon frequency is higher man the pump photon 
111, which separates the source signal into two iUumination 25 frequency) stimulated Brillouin scattering, stimulated Ray- 
source signals, one being transmitted to the reference probe le igh scattering, stimulated Rayleigh-wing scattering, four- 
108 and the other being transmitted to the sample probe 110. wavc nixing, and others which are well known to those 
The reference probe 108 focuses its iUumination source practiced in the art In typical stimulated scattering 
signal to the mirror 70* mounted on the back of the reference experiments, the pump photons are provided by an intense 
mirror 70 which is mounted on the galvanometer, and the 30 laser pump source, and the probe photons are either provided 
sample probe 110 transmits its illumination source signal to bv a weak probe beam or are provided by incoherent 
a fixed mirror 112. The interferometer also includes a scattering processes or noise The operation of a laser, for 
photodetector 114 for receiving the combination of light example is based on stimulated emission of radiation at the 
reflected back from the reference mirror 70 and the fixed lascr oscillation frequency which builds up from optical 
mirror 112, and for producing an analog signal 115 corre- 35 noisc Present in the laser cavity due to a spontanous (Le, 
sponding to the intensity of ligjit received Because a long- incoherent) emission background In stimulated Raman scat- 
coherence length Mumination source 106 is used the analog tering experiments, the probe beam may either be supplied 
mterferometric signal 115 produced by the photodetector as a weak source at a frequency cccrespondingtomeRaman 
114 wiU be a relatively constant amplitude sinusoidal signal transition to be interrogated in the sample or it may also be 
having a frequency equal to the Doppler shift corresponding 40 ^wed to build up from incoherent spontaneous Raman 
to velocity fluctuations in the reference mirror 70 experi- scattering noise In either case the intense pump radiation 
enced by the electric field in the reference arm. sets up a condition in which the probe experiences 

The analog signal 115 produced by the photodetector 114 frequency-dependent coherent gain in the medium. The 

is sent to an interval-detect circuit 116, for detecting features frequency dependence of the gain is determined by the 

in the signal 115 mat are regular in time (such as zero 45 mc <liHm , s specific atomic and molecular composition; thus 

crossings). These features 118 are fed into a clock generator probing the frequency dependence of the stimulated gain 

circuit 120 for generating a digital clock source signal 122 rnay serve as a sensitive probe of tissue biochemistry with 

for clocking (triggering) the analog-to-digital converter applications in medical diagnostics. In typical experiments 

device 82 used in the deconvolution scheme 80. of ^ rv P c * the probe radiation is monochromatic and the 

Accordingly, the sampling rate of the analog-to-digital con- 50 S^in experienced in traversing an excited medium is mea- 

verter 82 will be synchronized according to the fluctuations sure d using either direct, gated or synchronous detection 

in the reference mirror 70 translation velocity detected by techniques as the frequency of the probe radiation is 

the calibration interferometer 48. Examples of interval- scanned The primary idea of stimulated coherent spectro- 

detect and clock generator circuits for use with the present scopic OCT is to take advantage of the broad spectral 

invention include Tektronix 465 oscilloscope, commercially 55 extent of OCT probe light in combination with a separate 

available from Tektronix, inc. pump Ught beam to perform stimulated coherent spectros- 

As will be apparent to those of ordinary skill in the art, °°Py ovcr the whole source spectrum at once, while simul- 

there are many way to use the calibration interferometer 48 taneousiy using the short coherence length of the OCT probe 

to compensate for the inevitable velocity fluctuations in the li 8 nt to depth resolve the resulting stimulated scattering 
reference mirror speed; all of which are within the scope of 60 spectrum. 

the present invention. For example an alternate method for Referring again to FIG. 5, in stimulated coherent spec- 
incorporating the narrow-band iUumination source calibra- troscopic OCT an intense pump laser 130 is directed to the 
tion interferometer into the data acquisition system includes sample 66 at the appropriate frequency to induce depth- and 
the steps of: digitizing both the calibration and the SLD frequency-dependent gain in the sample volume interro- 
tnterf erograms at a sampling rate mat is higher than twice 65 gated The depth resolved sample scattering spectrum C(k) 2 
flie frequency of the calibration interferogram; detecting or sample arm power spectrum S„(k) obtained according to 
regular features, corresponding to regular intervals of space, the steps discussed above wfll thus contain features corre- 
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sponding to the frequency-dependent round-trip gain expe- 
rienced by the OCT source radiation (inelastic backscatter- 
ing characteristics of the scatterers resident within the 
window). Specific means for combining the pump and probe 
beams in SC-SPOCT are illustrated in FIGS. \\a and \\b. 

For stimulated scattering processes which do not require 
phase matching between the pump and probe beams, the 
pump beam may be directed into the sample at any angle 
with respect to the probe beam A convenient design is thus 
to combine the pump and probe beams coaxiaUy, as illus- 
trated in FIGS. \la and Ub. In the bulk optic SC-SPOCT 
system illustrated in FIG. 11a, the pump light 130 is com- 
bined with the sample arm light 204 in beamsplitter #2 206 
which may be either a broadband, wavelength-independent 
beamsplitter (similar to beamsplitter #1). or else preferably 
may be a dichroic beamsplitter which preferentially reflects 
the pump laser light and passes the OCT probe light It 
should be obvious that beamsplitter #2 206 may be placed at 
any point in the sample arm of the interferometer either 
before or after the sample arm focusing optics. Also, beam- 
splitter #2 may also be placed before beamsplitter #1 54, 
although this is not a preferred embodiment since interfer- 
ence between the long-coherence-length pump light in the 
sample and reference arms will make detection of the probe 
light interference more difficult Several other modifications 
to the normal OCT setup are required for SC-SPOCT. Since 
some of the intense pump light may be reflected from the 
sample and hence re-coupled into the interferometer, some 
optical elements designed to protect sensitive interferometer 
components and allow for optimal performance may be 
added. These may include an optical isolator 212 placed 
immediately after the OCT broadband source, and a 
wavelength-selective filter 216 placed before the detector 
74 . The wavelength-selective filter 216 may be a wavelength 
long-pass filter in the typical case in which the pump laser 
frequency is higher than the probe light frequency, or else it 
may be a wavelength short-pass filter in the opposite case or 
a interference band-pass (at the OCT probe wavelength) or 
band-reject (at the pump laser wavelength) filter in either 
case. It should also be noted mat the sample arm optics 
depicted in FIGS. 11a and 116 are only a schematic repre- 
sentation of whatever sophisticated optics may be present in 
the actual interferometer, including, for example, the relay 
coupling and focusing optics of a biomicroscopy 
endoscope, catheter, or other medical diagnostic device. 
Finally, a timing trigger 220 may be provided for the 
purpose of gated or synchronous detection (described later). 
This trigger may be derived from an independent trigger 
source, or from the trigger electronics of either the pump 
laser or OCT probe source. The trigger may also be obtained 
from detection of a small pick-off beam from either the 
pump laser or OCT probe light source, in a manner which is 
familiar to those skilled in the art 

FIG. lib depicts the coaxial combination of the pump and 
probe light in the case of a fiber optic OCT system, in which 
all of the extra elements are present in their fiber optic 
implementation but have the same role as their counterparts 
in the bulk optic system of FIG. Ha. In this case, the fiber 
optic version of the dichroic implementation of beamsplitter 
#2 is know as a wavelength division multiplexer (WDM). 
FIG. 11c illustrates the case in which the pump laser light 
and the OCT probe light are not delivered to the sample 
coaxialiy. but rather are delivered through either the same or 
separate optical focusing systems in such a way that they 
arrive with an angle 0 between them This arrangement may 
be used when the coaxial system of FIGS, Ma and lib is not 
practical, or else when a specific angle between the pump 
and probe beams is required to meet a phase matching 
condition. 
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The concept and process of SC-SPOCT is illustrated 
schematically in FIG. \2a~\le . OCT probe light 50 with a 

source spectrum given by S„ (k) is incident in the source arm 
of the SC-SPOCT interferometer. This light may have an 
5 arbitrary spectral shape, but it is characterized by a center 
spatial frequency ko and a full- width-half-maximum 
(FWHM) spectral bandwidth Ako. Light from the pump 
source 130 which is incident on the second beam splitter is 
also characterized by its center frequency kj and FWHM 
10 spectral bandwidth Ak^ In the preferred embodiment of 
SC-SPOCT. Ak t will be much narrower than Ak^ and in fact 
Ak t will set the maximum spectral resolution with which the 
stimulated coherent scattering features of the sample will be 
resolved. The spectrum of light incident on the sample will 
15 thus be a superposition of the source and pump spectra, as 
illustrated in the figure. The sample is characterized by its 
own intrinsic stimulated scattering spectrum C(k) 2 which is 
a function of the pump wavelength k x and the internal 
atomic/molecular makeup of the sample. In the example of 
20 FIG. lid, the stimulated scattering spectrum is characterized 
by several separate coherent scattering peaks corresponding 
to internal electronic or vibrational/rotafional transitions of 
the atoms and molecules within the sample, each peak 
having its own center frequency and spectral bandwidth. In 
25 the case of stimulated emission SPOCT, the stimulated 
scattering spectrum will reflect electronic energy states 
comprising the fluorescence spectrum of the sample having 
been stimulated with light at the spatial frequency k^ In this 
case, the stimulated scattering peaks will be quite broad 
30 (typically several hundred nanometers bandwidth in the 
wavelength units typically used in atomic spectroscopy), in 
many cases as broad or broader than the spectral bandwidth 
of the OCT probe source itself. In the case of stimulated 
Raman scattering, the stimulated scattering spectrum will 
35 reflect the Raman spectrum of the sarnie as a function of the 
frequency offset from the excitation frequency k x . and the 
bandwidth of the peaks will typically be on the order of 
single cm -1 to tens of cm -1 , in the spatial frequency units 
typically used in molecular spectroscopy. 
40 As illustrated in FIGS. Ma-Ale . the pump laser beam at 
frequency k x will induce a gain with frequency dependence 
given by C(k) 2 in the sample, and thus the light collected by 
the sample arm optics the sample will contain a combination 
of dastkally and inelastically scattered light whose com- 
45 posite spectrum will appear as illustrated. This composite 
spectrum will have the overall shape of the source spectra m, 
however for those frequencies in the stimulated spectrum of 
the sample C(k) 2 which overlap the source spectrum, stimu- 
lated gain will lead to alterations in the collected spectrum 
50 reflecting the frequency-dependent gain of the sample. The 
stimulated spectrum of the sample C(k) 2 and the power 
spectrum of the light returning from the sample S„(k) may 
both be recovered as a function of depth in the sample using 
the methods described in the previous sections of mis 
disclosure. 

In general, the spectrum of light returning from the 
sample S„(k) and the sample backscatter spectrum estimate 
C(k) 2 will include alterations due to the elastic basckscatter 
spectrum of the sample as well as due to stimulated scat- 
tering processes. Methods for isolating the stimulated scat- 
tering spectrum with high sensitivity by using gated or 
synchronous detection are described in the next paragraph. 
Assuming that the alterations to the backscatter spectrum are 
either small compared to alterations due to stimulated coher- 
ent scattering, or else have been otherwise removed, a 
method for extracting a quantitave parameter for the stimu- 
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iated scattering spectrum from estimates of S ss (k) is next interferometer. In the computer 234, these digitized values 
described. Here we extend the notation for SJk) to S„(l.k) m processed (using the methods described above) to gen- 
to denote the sample arm power spectrum estimated at a erate estimates of C(k) 2 and/or S w (Lk) both with and without 
distance into the sample of 1. We denote the frequency- the pump laser on. The estimate of C(Lk) and/or S J4 (l.k) 
dependent gain of the sample at a depth 1 as G<lic). and note 5 obtained with me iaser on is subtracted from the estimate 
that for stimulated coherent scattering processes the sample obtained with the laser off, generating a high sensitivity 
arm light will experience exponential gain between any two measurement of C(k) 2 and/or S„(l±) due only to coherent 
depths I, and 1 2 in the sample as a function of Gil Jrt and the scatteiin g Processes which is insensitive to the elastic back- 
round-trip propagation distance I 2 -i x : |ft sc f> T p^V,. * i u 

2 1 io Synchronous detection of stimulated coherent scattered 

^d^A^^Ga^WXi,-!,)] (AA) U 3 h \ is m nG ' lSa ~*- synchronous detection. 

the intensity of the pump laser light is modulated at a 

Here Gil^X) denotes the stimulated scattering gain speci- Modulation <* * e P um P Iff «y be 

fied over the depth range from 1, to L accomplished either sinusoidaUy by modulating the pump 

FIG. 13 illustrates how the ScWt sample arm power 15 ^ *»*[• * ° f * 1 ectow ^ c ° f 

c^^ n ,;iik. a'k^a * ^ ^ ^ ^ acousto-optic modulation. Alternatively, the pump laser 

TT^Zt ? ? ° f P r0 P^ ati0n *** light ma/ be modulated with a square wave^litude 

in a sample with stimulated scattering. Estimates of profile by use of an optkal c ^ ^ synZLs 

GKliU^Jc) may be obtained from estimates of S^(Lk) by detection, the detector 74 signal and a trigger 220 signal with 

solving equation AA: ^ a fixed phase relationship to the pump laser modulation are 

provided as signal and reference inputs, respectively, to a 

(hi i 1 lJ 5 ^ 2 ' k) ] (BB) lodc_m amplifier 236. The lock-in amplifier 236 outputs a 

2 ' ik-h)[S^ k )\ F°P°*tional to the content of the detector signal 

which is modulated at the same frequency as and is in phase 

25 with the pump laser modulation. The lock-in amplifier 

The actual stimulated coherent gain experienced in tissues output is digitized one or more times per coherence length 

in most cases will be very small, perhaps altering the sample of travel of the reference arm, and estimates of C(k) 2 and/or 

"™ £75 ^ as lMe as 1 m 10 3 or even 1 S s fiX) both with and without the pump laser on are obtained 

part in 10 . In addition, as mentioned above, the stimulated as for gated detection 

scattering spectral alterations will occur in addition to spec- 30 DC Stimnteted-Emission Spectroscopic Optical Coherence 

tral alterations due to elastic backscatter. It would thus be Tomography (SE-SPOCT) 

very desirable to have a method to isolate very small Stimulated emission spectroscopic optical coherence 

^SSS^t^ ^ A . tomography (SESPOCT) is a specific implementation of 

Methods for using gated or synchronous detection to stimulated coherent spectroscopic OCT which obtains a 

ac^eve this objective are fllustoted in FIGS. 14a-* and 35 depth-resolved stimulated emission spectrum of a sample. 

^/P e "nderiying idea of these methods is to employ SE-SPOCT may be used to image uptake of high quantum 

modulation of the pump light sourcepowerto modulate only efficiency laser dyes or infrared^mitting exogenous fluo- 

fcose component of the estimated sample arm power spec- rescent probes into biological specimens. Examples of the 

Hum which are due to stimulated scattering effects. latter include probes for measurement of mtrTcellular pH 

^ted detection of stnnulated coherent scattering is fllus- 40 (eg. carbocyanine. with an excitation rnaximum at 780 nm 
trated schematicaUy in FIG. 14a. For gated detection, a and an emission maximum at 795 nm), and nucleic acids 
puked pump laser is used. The pump laser light may be (e.g. IR-13Z with its excitation maximum at 805 nm and its 
pulsed either by using a Q-switched or mode-locked laser. emission maximum at 835 nm). A suitable pump laser 130 
Gated detection is the optimal detection method for stimu- for use with such dyes includes a 780 nm nan^econd^pulsed 
lated coherent scattering processes for which the stimulated 45 fiber-coupled diode laser. When such a pump laser is 
gain coefficient G<Uc) is a function of the pump laser utilized, an 830 nm superiuminescent diode is used as the 
intensity, Le, for nonlinear optical processes such as stimu- OCT probe source 50. Another pump laser 130 suitable for 
lated Raman scattering. For higher-order nonlinear optical use with these or other infrared-emitting fluorescent species 
stimulated scattering processes, it may also be advantageous is a NdrYAG^umped dye laser. When such a pump laser is 
to use a pulsed OCT probe source synchronized to the pump so utilized, a modelocked Ti:Al 2 0 3 laser may be used as a very 
laser pulse frequency, such as the broadband output of a broadband OCT probe source 50. The latter combination of 
femtosecond mode-locked laser. In gated detection, the pump and probe sources forms a broadly tunable system 
pump laser is triggered by a timing trigger source which may capable of detecting both depth- and frequency-dependent 
be either a separate system trigger, or may be a part of the gain in a range of clinically viable probes, 
pump laser pulse electronics. The detector 74 output is then 55 X. Raman Scattering Spectroscopic OCT(SRS-SPOCT) 
processed by a gated integrator device 224, which integrates Similar to SErSPOCT. the present invention also provides 
the detector output only for the duration of a gate which is a system or method for performing stimulated Raman scat- 
set to coincide with the laser pulse duration. The gated tcring spectroscopic OCT (SRS-SPOCT). In this 
integrator outputs a signal 226 representing the average of embodiment, the system or method includes the means for, 
the signal detected during the gate, optionally additionally 60 or the step of directing a high intensity pump light into the 
averaged over many separate laser pulses. The gated inte- sample interaction region. Resonant gain experienced by the 
grator also outputs a signal 228 proportional to the average low-coherence probe light will appear as localized peaks in 
of the signal detected during a separate gate signal when the the depth-resolved backscatter spectrum of the sample. The 
pump laser is not on. also optionally averaged over many SC-SPOCT concept is well suited to Raman spectroscopy 
separate such gates. These two outputs are then digitized in 65 since SRS spectra have sharp spectral features, and can be 
A/D converters 230. 232, one or more times per coherence obtained in any desired wavelength range by selection of the 
length of travel of the reference arm of the SC-SPOCT pump laser frequency. Thus substantial vibrational/ 



